16.451 Lecture 21 Smallest nucleus: The deuteron Nov 22, 2005 


Basic properties: J 

mass: me 2 = 1876.124 MeV 

binding energy: B = - M = m p + m n - m d = 2.2245731 MeV 

(measured via y-ray energy in n + p d + y) 
RMS radius: 1.963 ± 0.004 f m (from electron scattering) 

J n , I = 1 + , 0 (lectures 13, 14) 

[i - + 0.8573 |i N 


+ y) ^ 


quantum numbers: 
magnetic moment: 

electric quadrupole moment: Q 


+ 0.002859 ±0.00030 bn 

(-» the deuteron is not spherical! ....) 


Important because: 

• deuterium is the lightest nucleus and the only bound N-N state 

• testing ground for state-of-the art models of the N-N interaction. 


Electron scattering measurements: 



Q (fm' 1 ) 


Because the deuteron has spin 1, 
there are 3 form factors to describe 
elastic scattering: the "charge" (G c ), 
"electric quadrupole" ( G q ) and "magnetic 
(G m ) form factors. (JLab data) 


Combined Data 


Intrinsic Shape of the Deuteron 



Interesting feature: strong attractive np force, but 
a void in the center - the deuteron is hollow! ... Why? 
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Interpretation of quantum numbers: 


j*,i = r, o 
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S - S n + s p 
J = S + L - I 

k = (+) (+) (-1) L = + => I = 0, 2, 4 


• Of the possible quantum numbers, L = 0 has the lowest energy, so we expect 
the ground state to be L = 0, S = 1 (the deuteron has no excited states!) 

• The nonzero electric quadrupole moment suggests an admixture of L = 2 

(more later!) 


introduce Spectroscopic Notation: 2S+1 T 

^ J 

with naming convention: L = 0 is an S-state, L = 1 is a P-state, L = 2: D-state, etc... 

o 

the deuteron conf iguration is primarily 5^ 


Isospin and the N-N system: 
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i = 0,1 


The total wavef unction for two identical Fermions has to be antisymmetric w.r.to 
particle exchange: 

^ total ~ V space X 0spin x X isospin 


Central force problem: 

space = f(r ) Y lm {0,(/>) 

with symmetry (-1) L given by the spherical harmonic f unctions 


Spin and Isospin conf igurations: 

S = 0 and 1 = 0 are antisymmetric ; 
S = 1 and 1=1 are symmetric 


3 c state can only 
1 be I = 0 ! 
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Magnetic Moment: \i d = 0.857 pi N 


In general, the magnetic moment is a quantum-mechanical vector; it must be aligned 
along the “natural symmetry axis" of the system, given by the total angular momentum: 

jU ~ J 

But we don't know the direction of J , only its "length" and z-projection as expectation 
values: 

(j 2 ) = J(J + 1); (J z ) = mj = (-J.... + J) 


in a magnetic field, the energy depends on mj via A E =-(/}• B) = - Qj rrij B ju N 


Strategy: we will define the magnetic moment by its maximal projection on the 
z-axis, defined by the direction of the magnetic field, with mj = J 



Use expectation values of operators to calculate the result. 


Calculation of \i : 



M = mj=J = 9j J Mn 


Subtle point: we have to make two successive 
projections to evaluate the magnetic moment 
according to our definition, and the spin and 
orbital contributions enter with different weights. 


1. Project onto the direction of J: 


2. Project onto the z-axis with mj = J: 


/}• J = 

cos [3 = 


ll • J 

/d cos a = . 

jj(J + 1 ) 


mj_ = J_ 

i ji V J ( J+1 ) 


ju = g j J jj n = ji cos cc cos f3 


jU • J 


1 

(J+l) 


Next, we need to f igure out the operator for ju 
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spin and orbital contributions: 


We already know the intrinsic magnetic moments of the proton and neutron, so these 
must correspond to the spin contributions to the magnetic moment operator: 


Mp = 

+ 2.79 /% - g sp S 

Mn 

9 s,p 

+ 5.58 

Mn = 

-1.91/% = g sn S 

Mn 

9 s,n ~ 

-3.83 


For the orbital part, there is a contribution from the proton only, corresponding 
to a circulating current loop (semiclassical sketch, but the result is correct) 


M 


// = I n r £ = g £ £ jU N , 


For the deuteron, we want to use the magnetic moment operator: 


M ( 9 s,p ^ p 9s, n ^ n L p J Mn 



Details: 


1. L p — j L because m n = m p , and L is the total orbital angular momentum ! 

2. L = 0 in the "S-state" ( 3 Si ) but we will consider also a contribution from 

the “D-state" ( 3 Di) as an exercise 

3. The proton and neutron couple to 5 = 1, and the deuteron has J = 1 


M = \{ 

£- 3 ) = I 1 

[ (fls.p $p + g s ,n + ^l) • J ^ 

Mn 

Trick: use 

Sn + Sp — s 

and write the operator as: 




^ ^ ^ ^ 


05 

^ Ig 

II 

s,p 9 s,n ) ^ 

+ l(S s ,p^Ss!^(Sp-S n ) + \L j 

^ "" ^ 

) Mn 


But the proton and neutron spins are aligned, and • J / 

so the second term has to give zero! 


= (S n 
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Comparison to experiment: 


10 


Md = °- 857 Mn 


This is intermediate between the 
S-state and D-state values: 


li ( — 0.880 jUpj 

jU ( 3 D, ) = 0.310 jU N 


Suppose the wave function of the deuteron is a linear combination of S and D states: 


Wd) = a 


3 D 1 ) with a 2 + b 2 =1 


Then we can adjust the coeff icients to explain the magnetic moment: 


Md = 0-~b 2 ) Mi%) + b 2 mC - DO 


b 2 = 0.04, or a 4% D-state admixture accounts for the magnetic moment ! 
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Deuteron Summary so far 
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Quantum numbers: (J 71 , 1) = (1\ 0) favor a 3 5 x conf iguration (with S = 1, L = 0) as 
the lowest energy n-p bound state 

Magnetic moment: |u = 0.857 jli n is 2.6% smaller than for a pure 3 5 1 state and is 

consistent with a linear combination of L = 0 and L = 2 components: 


Vd 


> 


a 


%) + b 



with 


a 2 + b 2 =1 and 


b 2 = 0.04 


If this is correct, we can draw two conclusions about the N-N force: 


deuteron: 



\ 


1. The S = 1 configuration has lowest energy (i.e., observed deuteron 
quantum numbers) , so the N-N potential must be more attractive 
for total spin S = 1 than for S = 0. 

2. The lowest energy solution for a spherically symmetric potential is 
purely L = 0, and the L = 2 wave f unctions are orthogonal to L = 0 wave 
functions, so a physical deuteron state with mixed symmetry can only 
arise if the N-N potential is not exactly spherically symmetric! 


Also: the “hole" in the middle of the deuteron means that at very short distances 
the N-N potential must be repulsive - the neutron and proton do not overlap! 


Deuteron Quadrupole moment etc... 12 


Formalism for electrostatic moments: 

{see also D.J. Griffiths, ' Introduction to E/ectodynamics' Ch. 3, and 16.369) 

• Recall that an electric charge distribution of arbitrary shape can be described by an 
infinite series of multipole moments (spherical harmonic expansion), with terms of 
higher L ref lecting more complicated departures from spherical symmetry 

• The energy of such a distribution interacting with an external electric field E = - VV 
can be written as: 


dV 

1 d 2 V 


E mt = V(0) q + — 

Pz „ „ 2 

Qzz + ... 

oz 

o 4 3z 2 

0 


where the charge distribution is located at the origin and has a symmetry axis along z; 
q is the electric charge or ' monopole moment' 
p 2 is the electric dipole moment 
Q zz is the electric quadrupole moment 

and in general, higher order multipole moments couple to higher derivatives of the 
external potential at the origin. 
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Formalism, continued.... (axially symmetric case) 
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The multipole moments are defined in general by: E g 


with corresponding multipole operators proportional to 
spherical harmonic functions of order t 


= Jp(r)E, d 3 r 
E, = c t r‘ Y m (0) 


The f irst few multipole operators and associated moments are: 


order, i 

Moment 

Symbol 

operator, E £ 

0 

monopole 

R 

1 =V4 xr°Y 00 (6>) 

1 

dipole 

Pz 

z = rcosd = r 1 7 10 (^) 

2 

quadrupole 

Qzz 

r 2 (3cos 2 9 1 )= ^ r 2 Y 2O (0) 


Connection to nuclei and the deuteron problem: 
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For a nuclear system, the multipole moments are expectation values of the multipole 
operators, e.g. the electric charge: 

z 

q - j p(r) d 3 r - e £ j p'(r) E„ ^,(r) d 3 r = +Ze 

M / . 

The electric charge density is proportional to the probability density, i.e. 
the wave function squared, summed up for all the protons in the system !!! 


We already know the electric charge of the deuteron, but what about higher moments? 

integral 

Electric dipole moment, general case: must vanish! 



of the system space reflection t - 1 
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continued... 
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Basic symmetry property: no quantum system can have an electric dipole moment if its 
wave function has definite parity, i.e. if M 2 is even. 

(Incidentally, the possibility of a nonzero electric dipole moment of the 
neutron is of great current interest: so far the measured upper limit 
is |p n | < 10' 25 e - cm ....) 

In fact, all odd multipole moments must vanish for the same reason 


First nontrivial case: Electric Quadrupole Moment Q zz (F&H sec. 18.1) 

Q - =|M?) r 2 (3cos 2 «-l)dV 


ellipsoidal shape 


Q > 0 if b > a: “prolate" 
Q < 0 if b < a: “oblate" 



How to relate this to a nucleus, e.g. the deuteron? 
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Assume that the charge distribution is an ellipsoid of revolution with symmetry axis 
along the total angular momentum vector J: 



As for the magnetic dipole moment, 
we specify the intrinsic electric 
quadrupole moment as the expectation 
value when J is maximally aligned with z: 


Qint 

= 

< £ => 

II 

£ 



"Quantum geometry": 


cos p 


mj_ = J 

\j\ jjv+i) 


If an electric field gradient is applied along the z-axis as shown, the observable 
energy will shift by an amount corresponding to the intrinsic quadrupole moment 
transformed to a coordinate system rotated through angle (3 to align with the z-axis 
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Result: 
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Let Q )ab be the electric quadrupole moment we measure for the ellipsoidal charge 
distribution with mj = J: 



How do we apply this to anything? 

1. A spherically symmetric state (L = 0) has Q int = 0 (e.g. deuteron S-state) 

2. Even a distorted state with J = -§• will not have an observable quadrupole moment 

3. J = 1, L = 2 is the smallest value of total angular momentum for which we can observe 
a nonzero quadrupole moment 

-> these are the quantum numbers for the deuteron D- state! 


Deuteron intrinsic quadrupole moment: 
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recall our model for the deuteron wave function: 



a 




result for the quadrupole moment: 


Q lnt = ^ |a*b| (r 2 


note cancellation here 


/ 


SD 


— b 2 (r 2 

20 \ / DD 


= +0.00286 ± 0.00003 bn 


A good model of the N-N interaction can fit both the magnetic moment and the 
quadrupole moment of the deuteron with the same values of a and b! + 
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Basic features of the N-N potential, via the deuteron, etc.: 
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1. Independent of the value of L, the state with intrinsic spins coupled to S = 1 
has lower energy 

-> this implies a term proportional to: 




-1/4, S= 1 


+ 3/4, S = 0 


2. The deuteron quadrupole moment implies a non-central component, i.e. the potential 
is not spherically symmetric. Since the symmetry axis for Q is along J, Q > 0 
means that the matter distribution is stretched out along the J - axis: 



-> This implies a “tensor" force, proportional to: 



, ^ (Si •r) (s 7 •r) - - 

-(3 V 1 % V 2 - Sj-Sj 


Q > 0, observed 


compare: magnetic dipole-dipole 
interaction, see Griff iths problem 6.20 


N-N interaction continued.... 
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3. There is also a spin-orbit term, as deduced from N-N scattering experiments with 
a polarized beam: 

Vs-O ~ (i-s) 


(This plays a very important role also in determining the correct order of energy 
levels in nuclear spectra - more later!) 


4. Finally, all contributions to the N-N interaction are based on a microscopic meson 
exchange mechanism, as explored in assignment 5: 



Where M is a n, p, co ... meson, etc. 

and each term has a spatial dependence 
of the form: 

-m r 

e 

N 2 V(r) = g x(spin function ) 

r 
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State of the art N-N interaction model: (not to scare anybody...) 
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PHYSICS REPORTS (Review Section of Physics Utters) 149. No. 1 (1987) 1-89, Nonh-HoHand, Amsterdam 


(89 page exposition of one of only ~3 state-of-the-art models of the N-N 
interaction worldwide - constantly refined and updated since first release.) 


THE BONN MESON-EXCHANGE MODEL FOR THE 
NUCLEON-NUCLEON INTERACTION* 

R. MACHLEIDT* 

Los Alamos National Laboratory. MS H 85 Q , Los Alamos, NM 87545 , £/, 5, A . 
and 

Department of Physics, University of California, Los Angeles. CA 90024, U.S.A, 


and 

K. HOLINDE ++ and Ch. ELSTER 

Instimt fur Tkmretische Kernphysik tier Universtiat Bonn, Nussatiee 14 - 16 , D -5300 Bonn, West Germany 


Brace yourself - it takes 2 pages just to write all the terms down!!! 
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pseudoscater mesons: 


r ) = ^2 m n=[( ^f) o-,^+ z(«v) s t f(i 

(they use "a" for spin) 


(tensor operator) 


(F-6) 


scalar mesons: 


KK. r) = «,{[l - \ (^)] y K0 + ^ I V*Y(m t r) + Y(m,r)V*] + iZfrrfL-s] 


(F-7) 


vector mesons: 


(spin-orbit interaction) 


V,(m v , r) = m,|[l + ^ [F 2 K(m v r) + y(m,r)7 2 ] 

+ ^ (^) y(m v r)€r, ■ <r. - |Z,(«i I r)/.-S - iZ(m v r)S l2 } 

+ ^ Y(m,r)+ Y(m v r)tr , ■ tr, 

l ATT 

-4Z,(»i,r)£'S — |Z(m,r)S 12 } 
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more of this... 
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with 


Y(x) = z~*/x , Z(x) = (mjmf{l + 3lx + 3tx z )Y(x) * ^ 

„ , , , , , , iw, , r „ («t ‘ r K®i ' r ) 

2jW = )(llx + 1/-T )*W . 5 12 = 3 - ■■ - o\ • a- 2 , 


“Yukawa f unctions" 
and derivatives... 


(F.8) 


and 



Wc use units such that A - c = 1 (he- 197,3286 MeV fm), The use of the form factor, eq. (33), at each 
vertex (with n a = 1) leads to the following extended expressions: 





V.UUr) + 




VM..vr), 


(F.9) 


where , = A. + e, A_ , = 3, - s, gM, 1 . g = 10 MeV is an appropriate ch oice 


The full NN potential is the sum of the contributions from six mesons: 


VW- 2 V„(r) 


Parameters and Predictions: 24 


R. Machleidt et ai, The Bonn meson-exchange model for the nucleon-nucleon interaction 


Table 14 

Meson and low-energy parameters (LEP) for the confi guration space one-boson -exchang e 


potential (OBEPR) 

deuteron properties: 

UJs J 

m a (MeV) A„ (GeV) 

LEP 

Theory 

a 14.9 

138.03 

L3 

(MeV) 

2.2246 




P D (%) 

4,81 

p 0.95; [6.13 

769 

1.3 

Qt (fa 1 ) 

0.274 




fti (l^) 

0.8524 

T| 3 

548,8 

1.5 

A % (hn |;J ) 

(T8S50 




D/S 

0.0260 

a) 20; [0.0] 

782.6 

1.5 

r 6 (fm) 

1.9691 




a, (fa) 

-23.751 

B 2.6713 

983 

2.0 

r t (fa) 

2.662 




a, (fm) 

5.423 

o 7.7823 a 

550 s 

2.0 

r, (fm) 

1.759 

Impressively good agreement for ~ 

10,000 experimental 


low energy scattering 

data points in assorted n- 

-p and p-p scattering 


parameters, etc: 


experiments plus deuteron observables: y 2 /d.f <~ 1.07 ! 
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What does the NN potential look like? 
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It looks different in different spectroscopic states of the 2N system! 



Only the deuteron is bound! Its quantum numbers have the deepest potential well. 



